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1. Introduction. 



Gaudin's model describes a completely integrable quantum spin chain. Originally 
[1] it was formulated as a spin model related to the Lie algebra sl2- Later it was 
realized, cf. [2], Sect. 13.2.2 and [3], that one can associate such a model to any 
semi-simple complex Lie algebra g and a solution of the corresponding classical Yang- 
Baxter equation [4, 5]. In this work we will focus on the models, corresponding to 
the rational solutions. 

Denote by V\ the finite-dimensional irreducible representation of g of dominant 
integral highest weight A. Let (A) := (Ai,... , Ajv) be a set of dominant integral 
weights of g. Consider the tensor product V(x) := (/)... (E) V\j^ and associate with 
each factor of this tensor product a complex number Zi. The hamiltonians of 
Gaudin's model are mutually commuting operators = ^i{zi^ . . . , z^)^ z = 1, . . . , A^, 
acting on the space V(x)- 

Denote by (•, •) the invariant scalar product on g, normalized as in [6]. Let {/„}, a = 
1, . . . , J = dimg, be a basis of g and {/"} be the dual basis. For any element A G 
denote by A^'^ the operator l(x)...(x)A(x)...(x)l, which acts as A on the zth factor 



One of the main problems in Gaudin's model is to find the eigenvectors and the 
eigenvalues of these operators. 

Bethe ansatz method is perhaps the most effective method for solving this problem 
for g = sl2 [1]- As general references on Bethe ansatz, cf. [7, 8, 2, 9]. Applied to 
Gaudin's model, this method consists of the following. There is an obvious eigenvector 
in V(^x)- the tensor product of the highest weight vectors of the Va/s. One constructs 
other eigenvectors by acting on this vector by certain elementary operators, depending 
on auxiliary parameters Wi^ . . . , Wm- The vectors obtained this way are called Bethe 
vectors. Such a vector is an eigenvector of the Gaudin hamiltonians, if a certain 



of V(x) and as the identity on all other factors. Then 
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system of algebraic equations, involving the parameters Zi^s and tw^'s, and the highest 
weights Aj's, is satisfied. These equations are called Bethe ansatz equations. 

One usually constructs eigenvectors in statistical models associated to a general 
simple Lie algebra g by choosing a sequence of embeddings of Lie algebras of lower 
rank into g and inductively solving diagonalization problems for these subalgebras 
[15, 16, 3]. This leads to a rather complicated combinatorial algorithm, which very 
much depends on the structure of g. 

Another method was recently proposed by Babujian and Flume [37]. They con- 
sider as analogues of Bethe vectors, the rational V(x)-valued functions, which enter 
Schechtman-Varchenko solutions [30, 31] of the Knizhnik-Zamolodchikov (KZ) equa- 
tion [27] associated to g, cf. also [36, 38]. Varchenko and one of us show elsewhere 
[39] that such a vector can be obtained as a quasi-classical asymptotics of a solution 
of the KZ equation. This establishes a remarkable connection between a model of 
statistical mechanics and correlation functions of a conformal field theory. 

In this work we propose a new method of diagonalization of Gaudin's hamiltonians, 
which is based on Wakimoto modules over the affine algebra g at the critical level 
and the concept of invariant functionals (correlation functions) on tensor products 
of g-modules. This will allow us to treat the diagonalization problem and the KZ 
equation on equal footing and thus explain the connection between them. 

The critical level is the level k = —h^, where is the dual Coxeter number of 
g. The peculiarity of this value of level is that it is only for k = —K^ that the local 
completion f/fc(g)ioc of the universal enveloping algebra of g contains central elements. 
The center of U-h\i{g)ioc was described in [10, 11] (cf. also [12, 13, 14]). Elements of 
the center act on the space of functionals on tensor products of representations of g 
at the critical level, which are invariant with respect to a certain Lie algebra. One 
can interpret the hamiltonian as the action of a quadratic central element from 
f/_/iV (g)ioc on the zth factor of invariant functional. Other central elements give rise 
to other operators, which commute with each other and with Gaudin's hamiltonians. 
We call them higher Gaudin's hamiltonians. 

Wakimoto modules [17, 18] at the critical level are special bosonic representations 
of g, which are essentially parametrized by functions on the circle with values in 
the dual space of the Cartan subalgebra [) of g. We will construct eigenvectors of 
Gaudin's hamiltonians by restricting certain invariant functionals on tensor products 
of Wakimoto modules. In conformal field theory language, those are certain bosonic 
correlation functions. Analogues of Bethe ansatz equations will naturally appear 
as Kac-Kazhdan type [19] equations on the existence of certain singular vectors in 
Wakimoto modules. The existence of these singular vectors will ensure that the 
vector constructed this way is an eigenvector of the Gaudin hamiltonians and their 
generalizations. 

We will express the eigenvalues of these hamiltonians in terms of the Miura trans- 
formation and show that Bethe ansatz equations appear as certain analytic conditions 
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on the eigenvalues. The formulas we obtain for the eigenvalues generalize Sklyanin's 
formula [20] for g = sl2- The appearance of the Miura transformation in this con- 
text is not surprising, because the center of f/_/iV (g)ioc is isomorphic to the classical 
W-algebra associated to g^, the Langlands dual Lie algebra to g [10]. 

The space of invariant functionals makes sense for an arbitrary level k. For non- 
critical values of k it coincides with the space of genus correlation functions (or con- 
formal blocks) of the corresponding Wess-Zumino-Novikov-Witten (WZNW) model 
[21, 22, 23, 24, 25, 26]. This definition is equivalent to the more common definition of 
correlation functions as matrix elements of certain vertex operators, acting between 
representations of g. It is known that for non-critical k these functions satisfy the 
KZ equation [27, 28, 29]. 

One can obtain solutions of the KZ equation by restricting certain invariant func- 
tionals on tensor products of Wakimoto modules, but those of non-critical level. The 
rational function, which enters such a solution [30, 31] (cf. also [32, 33, 34, 26]), then 
appears as a bosonic correlation function, cf. [35], which coincides with the formula 
for an eigevector, obtained from Wakimoto modules at the critical level. This gives 
an explanation of the connection between the eigenvectors of the Gaudin model and 
solutions of the KZ equation associated to g, which was observed in [36, 37, 38, 39]. 

Gaudin's hamiltonians can be obtained from the transfer-matrix of a quantum 
spin chain related to affine quantum groups Uq[g) in a certain limit, when q ^ 1. 
A generalization of the results presented here to the case of Uq[g) will be published 
separately. 

We also want to remark that our construction of eigenvectors naturally fits into a 
program of geometric Langlands correspondence proposed by Drinfeld. This corre- 
spondence relates equivalence classes of G^-bundles over a complex algebraic curve 
£ with fiat connections, and D-modules on the moduli space of G-bundles over £. 

The paper is organized as follows. In Sect. 2 we recall the definition of Gaudin's 
model and Bethe ansatz procedure in the case g = sl2- Sect. 3 contains an inter- 
pretation of the model in terms of invariant functionals at the critical level. We 
show how singular vectors of imaginary weight give rise to Gaudin's hamiltonians 
and their generalizations. In Sect. 4 we recall the definition of Wakimoto modules 
and discuss the structure of the center of f/_/iV (g)ioc- We use these results in Sect. 
5 to construct eigenvectors of Gaudin's hamiltonians and to compute their eigenval- 
ues. In Sect. 6 we derive Schechtman-Varchenko solutions of the KZ equation, using 
Wakimoto modules. We show that the rational function entering a solution coincides 
with the formula for an eigenvector. In the Appendix we prove some technical results 
concerning Wakimoto modules, which are used in the main text. 

2. Gaudin's model. 

Let g be a simple Lie algebra over C of rank / and dimension J, and let U (g) 
be its universal enveloping algebra. For a dominant integral weight A denote by V\ 
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the irreducible representation of g of highest weight A. Denote by (•, •) the invariant 
scalar product on g, normalized as in [6]. Let /„, a = 1, . . . , J, be a root basis of 
g and be the dual basis. Denote by A the quadratic Casimir operator from the 
center of U (g): 

a = l 

The operator A acts on Va by multiplication by a number, which we will denote by 

Now fix a positive integer N and a set (A) := (Ai, . . . , Ajv) of dominant highest 
weights. Denote by V(x) the tensor product ® . . . ® ^Aj^- 

Let M, . . . , 2;jv be a set of distinct complex numbers. Introduce linear operator 
S'(m), which we will call the Gaudin hamiltonian, on V(x) by the formula 

(2.1) ^W = E-^ + E' ^ 



=1 - " 



where is given by formula (1.1). 

One can check directly [1, 3] that the operators Si,... , Sjv commute with each 
other. Therefore the operators S'(m) commute with each other for different values of 
u. (In the next section we will give a new proof of this fact.) 

One of the main problems in Gaudin's model is to diagonalize operators S'(m). This 
is equivalent to simultaneous diagonalization of the operators Si, . . . , Sjv- 

For g = sl2 this problem can be solved by algebraic Bethe ansatz [1, 2, 20]. 

Denote by v\ the highest weight vector of the module V\. The tensor product of 
the highest weight vectors 

|0) ■=vx,0...0 vxj, G V(x) 

is an eigenvector of the operators S^. Indeed, if and are not elements of the 
Cartan subalgebra of g, then the operator acts by on |0). If they are, then 

this operator multiplies |0) by a certain number. 
Therefore |0) is an eigenvector of S[u). 

The idea of the Bethe ansatz method is to produce new eigenvectors by applying 
certain elementary operators to the "vacuum" |0). 

Let {E^ be the standard basis of 5i2- Introduce operators F[w) on the space 

V(x) by the formula 

N p{i) 

(2.2) F{w) = ^ 



. ^ W — Zi 

8 = 1 ' 



where F^'^ is the operator, which acts as the generator F of s[2 on the zth factor of 
V(x) and as the identity on all other factors. Here w is a complex number, which is 
not equal to 2;i, . . . ^ z^. 
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Now consider Bethe 's vector 

(2.3) = F{wi) . . . F{w^)\0) 
in V(A). 

Explicit computation [1] shows that 

S{u)\wi, . . . ,w^) = s^{u)\wi, . . . ,w^) + 

(2.4) ^ I . 

j = l 3 

where Sra{u) is a function in u, and 

^ \- 9 

If all /j's vanish, then |t/;i,... ^Wm) is an eigenvector of S[u). The corresponding 
equations 

^ A- 2 

(2.5) ^ ^ Yl = 0' J = l,...,m, 

~; W,; — Zi ^ W,; — W. 

8 = 1 J ' Sfj J * 

are called Bethe ansatz equations. 

The eigenvalue of S[u) on the vector \wi^ . . . , 1/;^) is given by 

(2.6) s^{u) = ^Xmiuf - ^duXmiu), 
where 

(2.7) X™(«) = E-^-E 



■ U — Zi U — W;. 

8=1 ' J=l J 

There are several approaches to constructing eigenvectors of Gaudin's hamiltonians 
for general g. 

One can construct eigenvectors by choosing a sequence of embeddings of Lie al- 
gebras of lower rank into g and inductively solving the diagonalization problems 
for these subalgebras [15, 16, 3]. This leads to a rather complicated combinatorial 
algorithm, which very much depends on the type of g. 

Another method was proposed by Babujian and Flume in [37]. The eigenvectors 
should be constructed by applying to the vacuum |0), the operators 

N 

F,{w) = Y.^—. J = l,... 
~; w — Zi 

8 = 1 ' 

However, since such operators no longer commute with each other, one should also 
add some extra terms, which account for the commutators. The right formula can 
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be extracted from solutions of the KZ equation [30, 31] (and in fact can be obtained 
as quasi-classical asymptotics of such solutions [39]): 

(2.8) ..,<-)= Ht — ^|o)• 

Here the summation is taken over all ordered partitions U U . . . U of the set 
{1, . . . ,m}, where P = {il^i^^ ■ ■ ■ ^^ij}- Note that one can consider vector (2.8) as 
an element of the tensor product of Verma modules (/)... (E) with arbitrary 
highest weights Ai, . . . , Ajv- 

It was claimed in [37] that the vector \wi ^ . . . ,1^^) satisfies an analogue of the 
equation (2.4). Namely, 

(2.9) S{u)\w\\...,w'-) = s,,,...,Ju)H\ 
where ^^^^...^^^(m) is a rational function in u, 

■ -, W. — Zj ~r^- W;. — w. ' 

and Xj is a certain vector in V(x)- Thus, if the equations /^^ = are satisfied for all 
j = 1, . . . , m, then vector \w\^ , • • • , ^l^) is an eigenvector of Gaudin's hamiltonians. 
It is natural to call these equations Bethe equations and the vector \w\^ ^ . . . , tw^) a 
Bethe vector. 

The direct proof of formula (2.9) seems rather complicated, and it is not clear 
how to generalize it to higher Gaudin's hamiltonians, which we introduce in the next 
section. 

In this work we propose an alternative approach to the diagonalization of Gaudin's 
hamiltonians, which is based on a concept of invariant functionals (correlation func- 
tions) on tensor products of representations of affine algebras at the critical level. 

In the next section we will interpret Gaudin's model in terms of such functionals. 
We will show that Gaudin's hamiltonians are members of a family of commuting 
linear operators on the invariant functionals, which come from singular vectors of 
imaginary weight in the vacuum representation of g at the critical level. 

In Sect. 5 we will construct eigenvectors of Gaudin's hamiltonians by restricting 
certain invariant functionals on tensor products of Wakimoto modules. This con- 
struction will allow us to prove that if Bethe equations are satisfied, then vector 
\w\^ ^ . . . ,1^^) is an eigenvector of the operators S[u) and their generalizations, de- 
fined in the next section. We will also compute the eigenvalues of these operators. 
This construction will then be used in Sect. 6 to give an explanation of the appearance 
of Bethe vectors in solutions of the KZ equation. 



J=l J 



(2.10) 



f? = y: 
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3. Correlation functions. 

Let g be the affine algebra, corresponding to g. It is the extension of the Lie 
algebra g (/) by one- dimensional center CA^. Denote by the Lie subalgebra 

g(x)C[[t]]®CA^ of g. For any finite-dimensional representation Vx of g denote by the 
representation of on which g (/) tC[[t]] acts trivially and K acts by multiplication 
by A; G C. Denote by Y\ the induced representation of g of level k: 

Consider the projective line CP^ with a global coordinate t and N distinct finite 
points Zi^. . . , 2;jv G CP^. In the neighbourhood of each point Zi we have the local 
coordinate t — Zi. Denote g{zi) = g (/) — Zi)). Let gjv be the extension of the 
Lie algebra (SiLiSi^i) by one-dimensional center CK , such that its restriction to any 
summand g{zi) coincides with the standard extension. For f = [fi[t — Zi)^ . . . , fN{t — 
zn)) and g = {gi{t — Zi), . . . ,gN{t — zn)) from ®f^]g{zi) the cocycle cujv, defining 
this central extension, is given by 

(3.1) cuA.(f,g) = A^.^ — Uh.dg,), 

-, Ztti Ji 

where J- stands for the integral over a small contour around the point Zi. 

The Lie algebra g^ naturally acts on the tensor product V^;^^ = (E) ■ ■ ■ ® in 
particular, K acts by multiplication by k. We will say that we have assigned modules 
V|^, . . . , Yl^ to the points zi,... ,zn. 

Let 02 := gzi,...,ZN t>6 the Lie algebra of g-valued regular functions on CP^\{2;i, . . . , 
zn} (i.e. rational functions on CP^, which may have poles only at the points 
Zi^. . . ,^;jv), which vanish at infinity. Clearly, such a function can be expanded into 
a Laurent power series in the local coordinate t — Zi at each point Zi. This gives an 
element of g{zi). Thus, we obtain an embedding of g^ into the direct sum (SiLiSizi)- 
The restriction of the central extension to the image of this embedding is trivial. 
Indeed, according to formula (3.1), this restriction is given by the sum of all residues 
of a certain one-form on CP^; hence it should vanish. Therefore we can lift the 
embedding g^ — > (SiLiSizi) to an embedding g^ — > gjv. 

Denote by H^^~^ the space of linear functionals on V^_^j, which are invariant with 
respect to the action of the Lie algebra g^. Such a functional is a linear map : 
V^_^j — > C, which satisfies 

(3.2) i,{g.y) = 0^ V(/G0„J/G V(X). 

By construction, we have a canonical embedding of a finite-dimensional represen- 
tation Vx into the module Y\: 



X eVx 



1 (g) X G Va, 
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which commutes with the action of g on both spaces (recall that g is embedded into 
g as the constant subalgebra). Thus we have an embedding of V(x) = ® • • • ® V\j^ 
into V^_^j. We will keep the same notation V(x) for the image of this embedding. 

The restriction of any invariant functional G to the subspace V(x) C V^_^j 

defines a linear functional on V(x)- Thus we obtain a map : H^^~^ — )• V^\y 

Lemma 1. The map : H^^~^ — > V^\~^ is an isomorphism. 

Proof. Denote g+{zi) := g (/) C[[t — Zi]] C g{zi). We have a direct sum decomposition 

®Ls(^0 = (®Ls+(^0) ®Sz. 

Indeed, for g[t) G ® denote by g'^{t) and g~{t) its regular and singular 

parts, respectively. Any element f = {fi{t — Zi), . . . , /jv(t — zn)) of (SiLiSi^i) can be 
uniquely represented as the sum of 

N 

8 = 1 

and 

- z,) - - zi), • • • , - ^iv) - (r)iv(^ - ^n)) G ®Ls+(^o, 

where {f~)i{t — Zi) denotes the regular part of the expansion of the function /~ at 

Zf 

The gjv^module V^_^j is induced from the module V(x) over ®f-ig+{zi) ® CA^, on 
which ®f=i{t — Zi) ■ gj^{zi) acts trivially and K acts by multiplication by k. Therefore, 
as a g^-module, V^;^^ is isomorphic to the free module generated by V(x)- Hence the 
space of gz-invariant functionals on V^_^j is isomorphic to V^x^ ■ □ 

Lemma 1 shows that an invariant functional is uniquely defined by its restriction 
to the subspace V(x) C ^^fx)- Thus, for any rj G there exists a gz-invariant 

functional fj := (e^)~^?7 : V^_^j — > C, such that the restriction of fj to V(x) C ^^fx) 
coincides with rj. 

Consider the Lie algebra g° = 0z ® S C gjv, where g is the constant diagonal 
subalgebra of g^. Note that g does not lie in g^, because by definition elements of 
02 must vanish at infinity. One can consider the space of g°-invariant functionals 
on V^_^j. This space is isomorphic to the space of g-invariants of V(x)- In conformal 
field theory a g°-invariant functional on V^_^j is called a correlation function (more 
precisely, conformal block), and the equation (3.2) is called Ward's identity. 

Remark 1. Sometimes it is convenient to consider instead of the space of gz-invariant 
functionals on V^_^j, the space of coinvariants of Y{x) with respect to g^. These spaces 
are dual to each other. For a detailed study of the functor of coinvariants, cf. [25]. □ 
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Let M he a g-module, and let M* be the restricted dual linear space to M. We 
can define a structure of g-module on M* as follows: 

(3.3) [g ■ f]{m) = f{t{g) • F), / g M*,^ G g, m G M, 

where i stands for the Cartan anti-involution on g. It is defined on the generators of 
g as follows: 

t{Ei) = F„ t{Fi) = E,, i{Hi) = H,. 

This gives a structure of g-module on Af*, which is called contragradient to M. It 
is known that if the module M belongs to the category O of g-modules, so does 
M* . In particular, the module contragradient to an irreducible representation V\ 
with highest weight A is isomorphic to itself. Hence as a contragradient module, V^-^-^ 
is isomorphic to V(x)- Therefore H^x'j — ^(A)- 

Let Vq be the representation of g, which corresponds to the one-dimensional trivial 
g-module Vq. We will call it the vacuum module. Denote by Vq the generating vector 
of Vq. We assign the vacuum module to a finite point u G CP^, which is different 
from Zi, . . . , Z]\f. Denote by H^^ the space of gz^„-invariant functionals on V^;^^ ® Vq 
with respect to the Lie algebra gz,u- Lemma 1 tells us that H^\q) — ^(a)- 

Let X be a vector in Vq. For any rj G ^(\) consider the corresponding gz^„-invariant 
functional -q G H^xo)- restriction rj{-,X) to the subspace V(x) (E) X G V(a,o) defines 
another linear functional, ??'(•), on V(x)- Thus, we obtain a linear operator depending 
on M, X[u) : — )• which sends rj to rj' . Since V^■)^^^ ~ ^(a)? we can consider 

X[u) as an operator acting on V(x)- 

For A G and m ^TL^ denote by A(m) the element A (x) G 0. Now introduce 
the following vector in Vq: 

(3.4) s = \Y^U-\)r{-\w 

a = l 

This vector defines a linear operator 

S{u) : V^x) - V(x). 

Proposition 1. Tht operator S[u) coincides with the Gaudin hamiltonian. 
Proof. For any A G consider the element 

The expansion of g at Zi is equal to 



A 1 
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Therefore the action of g on the zth factor of the tensor product of V^;^qj is given by 



d 



n — l oo 



n 



i)\d 



m=0 



U — Zi 



\m + l 



Thus from (3.2) we obtain for any Y G V^_^j 



f3.5l 



fiiY,Ai-n)X) 



d 



n-l 



[n — 1)! dv/' 



\i=l m=0 



[U - Zi 



im + 1 ' 



This identity allows to "swap" the operator A{—n) from the module, assigned to the 
point M, to the modules, assigned to the points Zi, . . . , Z]\j. 

Put n = 1 in (3.5). Since any element of g (x) (t — 2;i)C[[t — Zi]] acts trivially on the 
subspace V(x) of ^\x}-, we obtain for any to G ^(xy- 



f3.6l 



rj [to 



/ N 

^(-1)^) = ^ E 
\i=i 



A(0 . 



to 



U - Zi 



Applying this identity twice to the element S G Vq, we obtain: 

(^ N d ja(i) N d t{j) \ 

^(-,5) = d-EE-— EE-^--,-o^ 



« = 1 a = l ^ J = l (1 = 1 -J 



N d 



z(/„)0-) z(/»)« 



-EE _ 

j = l a = l J 1 = 1 a = l ^ 



EE 



by formula (3.3). Hence 



N d 



^W = oEE 



+ EE 



2 , = 1 a = l - ^0' ^<J a = l " ^0(^^ " ^j) 

But this coincides with formula (2.1), because 

1 1 / 1 1 \ 



[u — Zi)[u — Zj) Zi — Zj \U — Zi U — Zj 



□ 



We can now prove that the operators S[u) commute for different values of u. In 
fact, we will prove a more general result. 



Let us specialize to the value k 



-h^ , where K" 



IS 



the dual Coxeter number. 



This value of level is called critical. For simplicity in what follows we will omit the 
superscript A;, when k = —K^. Thus, we will write Vq instead of Vq'' , etc. 

This value is special, because it is only at the critical level that S G Vq is a singular 
vector of imaginary weight. This means that A • X = for any A G ® C[t]. 
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Proposition 2. Let Z\ and Z2 be singular vectors of imaginary weight in Yq. Then 
for any pair of complex numbers u and v the corresponding linear operators Zi[u) 
and Z2{v) on V(x) commute. In particular, for Gaudin's hamiltonians we have: 

[S{u),S{v)] = 0. 

Proof. Consider a generalization of the previous construction, with two auxihary 
points, u and v. We can assign to each of these points the module Vq. The space 
H(xflfl) is again isomorphic to V(x)- Therefore we can define an action of a pair of 
elements, Xi,X2 G Vq, on V(x). Namely, for any rj G V(x) we consider the correspond- 
ing gz^„^^,-invariant functional rj G -ff(A,o,o)- Its restriction rj[-^ Xi^ X2) to the subspace 
^(A) ® Xi (x) X2 C V(x,o,o) defines another linear functional, ??'(•), on V(x)- Thus, we 
obtain a linear operator depending on u and (Xi,X2)(m,u) : V^^^ — > ^(a)? which 
sends rj to rj'. 

We will show now that if Xi and X2 are singular vectors of imaginary weight, then 
(Xi,X2)(m,u) = Xi[u)X2{v) and (Xi,X2)(m,u) = X2{v)Xi[u). This will prove that 
[Xl(u),X2(^;)] = 0. 

Indeed, Xi can be written as a linear combination of monomials ni) . . . 

Am{—n-m)vo^ where Aj G Q^n-j > 0. We can apply Ward's identity (3.6) to each 
summand and inductively "swap" operators n^) to the modules assigned to the 

points Zi^ . . . , 2;jv and v. This amounts to acting on them by 

in, - 1)! hi^- ^iP ^ - 1)! 9u-^-' - ^)"+V ■ 

But the last term of such a sum vanishes if we apply it to X2, by definition of singular 
vector of imaginary weight. Therefore this only affects the modules at the points Zi^ 
just as if we had Vq at the point v instead of X2. After that we can "swap" X2. 
We obtain: 

^(cu,Xl,X2) = [Xl(u)X2(^;)•??]H. 

Hence (Xi, X2)(m, v) = Xi[u)X2{v). In the same way we can show that (Xi, X2)(m, v) = 
X2{v)Xi[u)^ and Proposition follows. □ 

A description of the space Z[g) of singular vectors of imaginary weight in Vq is 
known [11, 10]. It is related to the structure of the center of the local completion 
f/_/iv(g)ioc of the universal enveloping algebra of g at the critical level [11, 10]. Let 
us recall this description. 

Recall that as a module over g_ = g (x) t~^C[t~^], Vq is isomorphic to its universal 
enveloping algebra f/(g_). Introduce a Z-grading on and on Vq by putting 

deg A(m) = — m,degUo = 0. There is an operator of derivative d = L^i of degree 
1 on the module Vq and hence on f/(g_), such that [5, A(m)] = —mA[m — 1), and 
d-vo = 0. 

Denote by Ji, . . . , J; the exponents of g. 
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Proposition 3. The space Z[g) of singular vectors of imaginary weight of Vq co- 
incides with the polynomial algebra in d^Si^i = 1,... > (applied to vector 
Vq), where Si^i = 1,... ,/, are mutually commuting elements o/f/(g_) of degrees 
+ 1, z = 1, . . . , /. 

A more precise description of the center will be given in the next section. Note 
that is equal to S. 

Proposition 3 and Proposition 2 show that S'i(m), . . . , Si[u) constitute a family of 
mutually commuting linear operators on the space V(x)- We call Si[u)^i > 1, higher 
Gaudin's hamiltonians. It would be interesting to find explicit formulas for them. 

Remark 2. For Z G Z[g) denote = Res„=^^(M — Zi)'"^Z{u). These operators 

generalize the operators = given by (1.1). They all mutually commute. 

Alternatively these operators can be defined as follows. 

Since Vq is a vertex algebra [48], we can associate a local current X{z) to any 
X G Vq. Denote by the (— m — l)st Fourier component of this current. This 

is a central element of f/_/iv(g)ioc [10]. Therefore it acts on the zth component of the 
tensor product Vaj ... Va^^. The dual operator gives rise to an operator on the 
space of gz-invariant functionals, -ff(A) — ^(A)- Using the generalized Ward identity 
from Sect. 5, we can show that the corresponding operator on V(x) coincides with 

Note that we can define the space of gz-invariant functionals in a more general 
situation. To any g-module M we can associate the induced g-module of level A;, 

M = f/(s) ®uc,,) M. 

Fix N representations Afi,... , Af jv of g, and let Mi,... , Mjv be the induced g- 
modules. We can define the space of gz-invariant functionals on the tensor product 
Ml (g) . . . (g) Mjv. This space is isomorphic to (g) . . . (g) M^, where M* denotes 
the module contragradient to Mi. In the same way as above we can associate to a 
singular vector of imaginary weight Z G Vq, a family of commuting linear operators 
on this space. 

It is clear from construction that the operators Z[u) can be viewed as elements of 
the A^th tensor power U{g)^^ of the universal enveloping algebra of g, depending on 
Zi^. . . ,^;jv, and u. Therefore they automatically act on any A^-fold tensor product 
of g-modules. 

Denote by M\ the Verma module over with highest weight A G f)* over g, and by 
the contragradient module. 

From now on we will choose as our modules, M^^ with arbitrary G f)*,z = 
1, . . . , /. Then the space of gz-invariant functionals is isomorphic to the tensor prod- 
uct of Verma modules, ®^iAfA,. 
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If Aj's are integral dominant weights, then there is a surjectiveg® -homomorphism 
®i^i^A, ®i^i^A,- The operators Z{u), Z G 2{q), act on both spaces and commute 
with this homomorphism. Thus, the projection of any eigenvector of Z[u) from 
®i^i^A, on (E)iLiVx, is again an eigenvector. 

Our construction of eigenvectors of Gaudin's hamiltonians in the tensor product 
Ma J (x) . . . (x) will be based on Wakimoto modules. 

4. Wakimoto modules at the critical level. 

In this section we briefly describe important facts about Wakimoto modules [17, 
18, 40] (cf. also [41]). 

Let us first introduce the Heisenberg Lie algebra r(g). It has generators a„(m), 
a* (m), m G Z, a G A_|_, where A_|_ is the set of positive roots of g, and central element 
1. They satisfy the relations: 

[aa{n),a*p{m)] = [a„(n), a^(m)] = 0, [a* (n), aj^(m)] = 0. 

Now let M be the irreducible representation of r(g), which is generated from the 
vacuum vector satisfying 

(4.1) a„(m)u = 0,m > 0, a*(m)u = 0,m > 0, 

for all a G A_|_, and Iv = v. 

Let [) be the commutative Lie algebra [) (/) C{{z)). It has a linear basis hi[n) = 
Hi (/) t^^i = 1, . . . , /, n G Z. Any one-form x{^)dz G f)* ® C{{z))dz defines a one- 
dimensional representation c^(^) of [): 

^xW[/(^)] = ^ / Xiz)[fiz)]dz, fiz) G ^, 

where the integral is taken around the origin. 

Wakimoto modules constitute a family of representations of the affine algebra g of 
the critical level on the space M (/) cr^(^z)- 

In order to define these representations, we should construct a homomorphism p 
from g to the local completion of Ui[T[g)) (/) U{^)- Here Ui[T[g)) stands for the 
Heisenberg algebra, which is the quotient of the universal enveloping algebra of T[g) 
by the ideal generated by (1 — 1). We recall the construction of the homomorphism 
P- 

First we have to choose coordinates on the big cell U of the fiag manifold B-\G of 
g. This big cell is isomorphic to the nilpotent subgroup A^_|_ of the Lie group G of g. 
Since A^_|_ is isomorphic to its Lie algebra n+ via the exponential map, U is isomorphic 
to a linear space with coordinates G A_|_. We assign to the coordinate Xa the 

degree a. The right action of G on the fiag manifold gives an embedding of the Lie 
algebra g into the Lie algebra of vector fields on U. This embedding can be lifted to 
a family of embeddings into the Lie algebra of the first order differential operators. 
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depending on x ^ ^* ■ Denote these embeddings by p^. Exphcitly, we have for the 
generators Ei, Hi, Fi, z = 1, ...,/, of g: 

where Pp is a certain polynomial in Xa of degree /3 — a^; 

and 

Px{Fi)= E Qpj^ + x{Hi)x^,. 

where is a certain polynomial in Xa of degree /3 + a^. 

The homomorphism p-^ defines a structure of g-module on the space of algebraic 
functions on the big cell, C[a;„]„£A+- This module is isomorphic to the module Af*, 
which is contragradient to the Verma module over g of highest weight x- 

We also have another map from the Lie algebra n+ to the Lie algebra of vector 
fields on f/ ~ A^_|_, which comes from the Itjt infinitesimal action of n+ on its Lie 
group. The zth generator of n+ maps to the vector field 
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(^^Oi /JgAj. '^^(i 



where i?^ is a certain polynomial in Xa of degree /3 — a^. These operators generate 
another action of n+, which commutes with the action of n+, defined by p-^. 
Now we can construct the homomorphism p. Introduce the notation 

A{z) = Y.A{n)z-^-\ AGS 

and 

««(^;) = E 0'c,{n)z~'"'~^ , a^{z) = E ^al^)^"", M^) = E h,{n)z~'''~'^ . 

n^Wj n^Wj n^Wj 



We put 



p[E,{z)] = a„,(z) + E : Pp{z)«'p{z) 
Pmz)] = - E W ■■ o*f,{z)ap{z) : +h,{z), 
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where q is a constant. Here Pp{z) and Q^p{z) stand for the expressions, which are 
obtained from the polynomials Pp and by insering a*a{z) instead of and dots 
denote normal ordering. 

We also define operators Gi[n)^ z = 1, . . . , /, n G Z, by the generating functions 

(4.2) G,{z) = J2 G.{n)z-^-' = -a^,{z) + ^ : R}{z)ap{z) : . 

Theorem 1. There exist such Ci in the formula for p[Fi[z)]^i = 1,... , /, that the 
map p defines a homomorphism from the Lie algebra g to the local completion of 
f/i(r(g)) (x) U[[)), which sends the central element K G —h^. 

This was proved in [17] for g = s[2 and in [18, 40] in general. 

Example. We will write down the precise formulas in the case g = s[2 [17]. We will 
omit the unnecessary subscripts. We have: 

p[E{z)] = a{z), 

p[H{z)] = -2 : a{z)a*{z) : +h{z), 
p[F[z)] = — : a[z)a*[z)a*[z) : -\-h[z)a*[z) — 2dza*{z), 

and 

G{z) = —a{z). 
For explicit formulas in the case g = s[„, cf. [18]. □ 

For any x{^)dz G f)* ® C{{z))dz^ we can consider representation M (/) of 
r(g)) ® [). According to Theorem 1, the homomorphism p defines on M (/) a 
structure of g-module of the critical level, which we denote by W^(^). We call this 
module Wakimoto module with highest weight x{^)- 

Remark 3. The modules, which were originally constructed in [17, 18], only depended 
on X G f)* — [)*dz/z. However, the generalization of that construction to arbitrary 
x{z) is straightforward. □ 

Remark 4. The Laurent power series x{^) should be considered as an element of the 
space of — /i^-connections on a principal i7-bundle on the (formal) punctured disc, 
where H is a Lie group of [). Locally such a connection can be written as a first order 
differential operator —h^d^ + x{z)dz^ where x{^) is an element of [) (x) C{{z))^ which 
is isomorphic to [)* (/) C{{z))dz via the non-degenerate scalar product (•, •) on [)*. We 
can view such a connection as a linear functional on [), which is equal to —h^ on the 
central element. Here [) is the Heisenberg Lie algebra, which is the central extension 
of the Lie algebra of gauge transformations of this bundle. 

The space of connections is a torsor over the space of one-forms [)* (/) C{{z))dz. If 
we trivialize the bundle, we can identify this torsor with the space of one-forms. □ 
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The following result, which is proved in the Appendix, will be important in the 
next section. 

Lemma 2. Let fi[z) be a highest weight of the form 

oo 

^ n=0 

The vector Gi[ — l)v G W^^(^) is a singular vector of imaginary weight, if and only if 

We can also define another representation of g at the critical level. Let ttq — 
<C[hi[n)]^i = 1, . . . , /, n > 0, be the [)-module, on which hi[n)^n > 0, act by and 
hi[n)^ n < act as on the free module with one generator. Then Af (x)7ro is a g-module 
via the homomorphism p. We denote this module by Wq. 

One can check directly that the generating vector v = v (/) 1 of Wo is annihilated 
by A(m), A G 0, m > 0. Therefore the map p defines a homomorphism p : Vq — ^ Wo, 
which sends Vq G Vq to u G Wq. One can show that p is a embedding. Under this 
homomorphism all singular vectors of imaginary weight from Vq get mapped into the 
subspace ttq C Wq [10]. Thus, we obtain an injective map pz '■ 2(g) — )• ttq. 

In [10, 11] the image of this map was described. Let us recall this description. 

We can consider the algebra ttq as the algebra of differential polynomials. The 
action of derivative d on ttq is defined by its action on the generators: dhi[n) = 
—nhi[n — 1), and by the Leibnitz rule. The action of d on Z[g) and on ttq commutes 
with the map pz- Therefore the image of Z[g) is a differential subalgebra of ttq. 

On the other hand, for any g one can define a remarkable differential subalgebra 
of TTo = 7ro(g), which consists of "densities" of the corresponding classical W-algebra, 
defined by the Drinfeld-Sokolov reduction [42, 43]. This subalgebra can be charac- 
terized as the space of invariants of the nilpotent subalgebra n+ of g, which acts on 
TTo, cf. [44], Sect. 2.4. It can be considered as a classical limit of the vertex operator 
algebra of the quantum W-algebra. We denote this subalgebra by yV(g). 

Now let g^ be the Langlands dual Lie algebra of g. Recall that the Cartan matrix 
of g^ is the transpose of the Cartan matrix of g. Therefore there is an isomorphism 
between the Cartan subalgebras of g and g^, which preserves the scalar products. 
This induces an isomorphism between 7ro(g) and 7ro(g^). Therefore we can consider 
yV(g^) as a differential subalgebra of 7ro(g). 

Theorem 2. [10, 11] The image of Z[g) in Tro{g) under the homomorphism pz co- 
incides with yV(g^). 

Proposition 3 follows from this Theorem, because it is known that yV(g^) is a 
free commutative subalgebra of ttq generated by d^Pi^i = 1,... > 0, where 

deg Pi = di-\-l [43, 44]. The generator Pi is the image of Si G Z[g) from Proposition 3. 
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Example. If g = s[2, then ttq = C[h[n)]n<o- The image of Z[g) in ttq consists of 
polynomials in c)™P, m > 0, where 



This space coincides with yV(s[2), cf. [44], Sect. 2.1. □ 

The center 3(s) of the local completion f/_/iV (g)ioc of the universal enveloping alge- 
bra of g at the critical level was described in [10]. It consists of "local expressions" in 
^i, . . . ,5;, i.e. Fourier components of all differential polynomials in Si[z)^ . . . , Si[z)^ 
where Si[z) is the local current associated to Si in the vertex algebra Vq. There- 
fore 3(s) is isomorphic to the classical W-algebra, associated to by means of the 
Drinfeld-Sokolov reduction, cf. [10]. This W-algebra is the space of local functionals 
on a certain hamiltonian space H[g^) [43]. 

Now consider the space J-'q of all Fourier components of differential polynomials 
in hi[z)^ . . . , hi[z)^ cf. [44]. It is isomorphic to the space of local functionals on the 
space T[g) of — /^^-connections —h^d^ + x{^)dz. The map pz '■ 2(g) — )• ttq gives us 
an embedding p-^ : 3(s) — ^ ^o- The corresponding map of "spectra" T[g) — )• H[g^) 
is nothing but the generalized Miura transformation [43]. 

Flements of 3(s) act on Wakimoto modules by multiplication by constants. One 
can describe the map p-^ in terms of Wakimoto modules as follows: for X G 3(s), 
the value of P'^iX) at —h^d^ + x{^)dz is equal to the action of X on the Wakimoto 
module W^(^). This shows that the map p-^ is an affine analogue of the Harish- 
Chandra homomorphism. 

Recall that the Harish- Chandra homomorphism is a map from the center 3(s) 
of the universal enveloping algebra U (g) of g to the algebra of polynomials on the 
dual space [)* of the Cartan subalgebra [) of g. Any element of 3(s) acts on Verma 
modules by multiplication by a constant. Harish- Chandra homomorphism maps a 
central element X to the polynomial on [)*, whose value at x is equal to the action 
of X on M^. 

Example. If = s[2, then 3(s) consists of all local expressions in the components Sn 
of the current 




Siz) = -Y.:Liz)riz):=Y.S^z- 




Then we have according to the previous example, 

s{z) = \Hzf-\dMz). 
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Therefore the central element s„ acts on the module W^(^) by multiplication by 

^ / ^ A A ^ 2 ^ ' 

where we put x{z) = Emez x^")^;"""^ □ 

In the next section we will use Wakimoto modules and the description of the center 
to construct eigenvectors of the Gaudin hamiltonians and to compute their spectrum. 

5. BETHE VECTORS FROM WAKIMOTO MODULES. 

We will construct eigenvectors using invariant functionals on tensor products of 
Wakimoto modules. We will restrict ourselves with modules W^(^), for which x{^) 
has the form 

^ n=0 

In other word, we consider only the modules associated to — /i^-connections on the 
formal disc, which are regular or have regular singularity at the origin. We will call 
such highest weight x{^) regular. 

Slightly abusing notation we will write x instead of x^~^^ = Res^=o x{^)- 
Let us recall the geometric construction of the Heisenberg algebra F := r(g), cf., 
e.g., [45]. Consider the spaces = U ® C((t)) and T'^ = U ® C{{t))dt of functions 
and one-forms on the formal punctured disc with values in the linear space U with 
coordinates G A_|_. There is a natural non-degenerate pairing between them: 

(•,•) : ^° X ^ C, 

which sends {f{t),g{t)dt) to [2Tr J {f[t)^ g[t))dt. Here (•,•) denotes the scalar 
product on the linear space f/, with respect to which orthonormal. To 

the scalar product (•,•) we can associate in the standard way a central extension 
of the commutative Lie algebra J-^ ® J-"^. This is our Heisenberg Lie algebra F. 
The generators aa{m)^ a*^[m) correspond to Xa ® t'"^ G and Xa ® t'"^~^di G J^^ , 
respectively. 

Let set of distinct points on the projective line. Consider 

Wakimoto modules W^^(^) with regular highest weights = 1, . . . ,p. 

Let us choose a global coordinate t on CP^ and the corresponding local coordinates 
t — Xi, z = 1, . . . , p, at our points. Denote V{xi) = U (/) C((t — Xi)) (BU (/) C((t — Xi))dt. 
Let Tp be the central extension of the Lie algebra (B^-iT[xi) by a one-dimensional 
center, which coincides with our standard central extension on each of the summands. 
The Lie algebra Tp acts on Af®'' in a natural way, in particular, the central element 
1 acts as the identity. 

Consider the commutative Lie algebra J-'^ := ® J^^^ where J-^ (-^x) 
space of f/-valued regular functions (one-forms) on CP^\{a;i, . . . ,3;^}, which vanish 
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(have regular singularity) at infinity. We have an embedding of Ty^ into 
obtained by expanding a function and a one-form at a given point. The restriction 
of the central extension to the image of this embedding is trivial. Therefore we can 
lift it to an embedding Ty^ — > Pp. 

We also introduce the commutative Lie algebra [)x of [)-valued regular functions 
on CP^\{a;i, . . . , x^}, which vanish at infinity. We have an embedding [)x — ^ ® 
C((t — Xi)). Denote "xiz) = X\{'^)-, - ■ ■ ■, Xv{'^)i where i^(^z) G f)* ® 'C,{{z))dz and let cr^(^) 
be the one-dimensional representation of ®f=if) ® C((t — Xj)), on which [) (x) C((t — Xj)) 
acts according to its character x«(-2^)- Thus, [)x also acts on c^(^). 

The Lie algebra Tix := -^x ® {)x acts on 

Denote by J-xiz) the space of Tix^nvariant linear functionals on this tensor product. 
Proposition 4. Suppose that 

. °° 

X^{Z) = — + 2^Xi Z , 1 = 1,..., p. 

Z n 

U Xi{t ~ ^i) is the expansion of 

i=i ^ 

at the point Xi, for all i = 1,... ,p, then the space Jx(z) is one- dimensional. It is 
generated by a functional, whose value on the tensor product of the vacuum vectors 
of Wy^^l^z) is equal to 1. 
Otherwise Jx(z) = 0. 

Proof. The space Jx(z) is isomorphic to the tensor product of the dual space of the 
space of coinvariants of M®'^ with respect to the action of J^x, and the space of 
[)x-invariants of c^(^). 

Since M is free over the Lie subalgebra of L, generated by aa{m), a G A_|_, m < 0, 
and a* (m), a G A_|_, m < 0, we can use the same argument as in the proof of Lemma 1. 
We conclude that the space of coinvariants of Af®'' with respect to the action of J-'^ 
is one-dimensional and that the projection of the tensor product of the vacuum 
vectors of W^^(^) on this space is non-trivial. Therefore there exists an Tix^nvariant 
functional on (x)f_;^ W^^(^) , whose value on the tensor product of the vacuum vectors 
is equal to 1. 

Let f[t) be an element of [)x. This is a meromorphic [)-valued function on CP^, 
which may have singularities only at x^'s and which vanishes at infinity. It acts on 
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the one- dimensional module c^(^) by multiplication by 



lx.it-zi)[fit)]dt, 

■ -, ZTTl Ji 



where stands for the integral over a small contour around Xi. This sum vanishes 
for any /(t), if and only if Xi{t ~ ^i) is the expansion of 



X 



X^ 



i=l ^ 



at the point Xi for all z = 1, . . . Therefore the space of [)x-invariants of cr^(^) is 
non-zero, if and only if the condition of Proposition is satisfied. □ 

Now fix highest weights of g, Ai , . . . , Ajv, and a set of simple roots of g, , . . . , a^^. 
Consider the function 

^ A- ™ a- 

Denote by Xi[t — Zi) the expansions of \{t) at the points Zi^i = 1, . . . , A^, and by 
— Wj) the expansions of \{t) at the points w^^j = 1, . . . , m. We have: 

A.(^) = ^ + ... , ^^^{z) = -'^ + ^if^ + 

z z 

where 



^ A,: ^ a; 



— Zi ^ t/;,- — 



- 



We put p = -|- m, and denote = Zi^Xii^) = ^il-^^),^ = , A^; xjv+j = 

tWj, Xjv-i-j(-2) = = 1, . . . ,m. According to Proposition 4, the space J\(z),ijl(z) 

of 7iz,w^nvariant functionals on the tensor product ®f=-[W\^(z) ®f=\ ^Atj(^) one- 
dimensional. It is generated by a functional Tjv,m, whose value on the tensor product 
of the vacuum vectors is equal to 1. This is a linear map 

(5.2) TN,.m : ®tiWx,(.) ^ C. 

We will obtain eigenvectors of the operators Z{u), Z G Z[g)^ by restricting the func- 
tional Tjv,m to a certain subspace. 

According to Lemma 2, the vectors Gi^[ — l)v G = 1, . . . ,m, are singular 

of imaginary weight, if and only if the following system of equations is satisfied 

5.3 2^ --2^ ^ = 0, j = l,...,m. 

These are precisely the Bethe ansatz equations. 
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Denote by W^a(^) the subspace of a Wakimoto module VFa(^), which is generated from 
the vector v by the operators a* (0), a G A_|_. This space is stable under the action of 
the constant subalgebra g of g, and is isomorphic to the module contragradient 
to the Verma module over g with highest weight A = Res^=o ^{^)- Indeed, only the 
0th components a„(0) and a*(0) will contribute to the action of g on VFa(^). Our 
formulas show that if we restrict the homomorphism p to g and replace a* (0) by Xa 
and a„(0) by d/dxa^ we will obtain the formulas for the homomorphism which 
defines on C[a;„]„£A+ — ^\{z) ^ structure of g-module isomorphic to M^. 

The restriction of the homomorphism Tjv,m (5.2) to the subspace 

(5.4) (g>l, Wx,i,) (X) G,,{-l)v (X) . . . ® G,J-l)v 

defines a linear functional 

The g-module (E)i=iMy is graded by the weights of the Cartan subalgebra [) of g, 
so that all homogeneous components are finite-dimensional. We will show in the 
proof of Lemma 3 that the functional t/j[wi ^ . . . ,1^^) vanishes on all homogeneous 
components except the one of weight J2iLi ~ YlT=i ■ Therefore it corresponds to 



a vector (/){w\^ w^m ) in the tensor product of the Verma modules over g, ®f-iM\^ 

of weight EL^« -E7=i«., . 

There is a differential algebra homomorphism rjv,m from ttq to the algebra of ra- 
tional functions in u, having singularities at Zi,... ^ z^^Wi^ . . . ^w^m on which the 
derivative acts as du- In order to define it we have to specify the values of this 
homomorphism on the elements hs{ — l) of ttq. We put 



N 

8 = 1 ' 1 = 1 



(5.5) rN,m{hs{-t)) = hs{u) 

This implies that 

rNAKini) . . . hsJuM)) = n (_^^_i)!g^-.,-i ^- 



Theorem 3. If Bethe ansatz equations (5.3) are satisfied, then vector (f>[w\^ ^ . . . ,1^^) 
is an eigenvector of the operator Z[u) for any Z G Z[g) with the eigenvalue r^^mPziZ) . 

In order to prove the Theorem, we will need some properties of the vertex algebra 
Wo, cf. [10, 11]. Each vector A G Wo defines a formal power series A[z)^ whose 
coefficients are linear operators acting on representations of r(g) ® [). We call these 
formal power series local currents (note that they were called _^eMs in [10, 11]). They 
satisfy all axioms of vertex operator algebras [46, 47], except for the existence of the 
Virasoro element; that is why we call Wo a vertex algebra. 
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The local current, corresponding to a monomial basis element Ai(ni) . . . Am{'nm)v G 
/q, where Ai stands for a„, a*, or hi^ is given by 



Here we put c/(a„) = = 1, J(a*) = 0. This current does not depend on the 

order in which we apply normal ordering. In general one should use normal ordering 
nested from right to left. 

The space Vq is also a vertex algebra (but not a vertex operator algebra) [48], and 
the homomorphism p defines an embedding of this vertex algebra into Wo . 

Let us assign representations Li of T[g) ® [) to the points x^, z = 1, . . . ,p, of CP^. 
Denote by J{Li^ . . . , Lp) the space of Tix^nvariant functionals on (E)^-iLi. 

Introduce a Z-grading d on Wo by putting J(a„(n)) = J(a*(n)) = d[hi[n)) = —n. 
Let Y[z) be the local current corresponding to an element Y G Wo of degree d[Y). We 
will identify the (1 - J(y))-differential {t - dt^-"^^^^ G C{{t - Xij)dt^-'^^^^ 

with the operator Y[nY^'^ acting on the module (E)^-iLi. Consider the space J^^-^(^) 
of regular (1 — J(y))-differentials on CP^\{a;i, . . . , x^}, which have zero of order at 
least 2d(Y) — 1 at infinity. Any element y{t) of this space can be expanded in powers 
of t — Xi at each point Xi. The corresponding formal power series y^''\t — Xi) can then 
be viewed as an infinite sum of operators J2m>-d{Y) Cm^(?7i)^'-' on the space (E)'i=iLi. 
Denote 

y = j2y^'Ht-xi). 

8 = 1 

Proposition 5. For any r G J{Li, . . . ,Lp),Y G Wo, and X G ®Li-^8; T{y-X) = 0. 

In the case when Y[z) is one of the elementary currents: a„(2;), a*^{z)^ or hi[z)^ this 
follows from the definition of the space of invariant functionals. The general result 
follows by induction on the power of Y. By linearity, it is sufficient to prove the 
statement in the case when y{t) = (t — Xj)""'"''^^)"^ n < —d[Y). Proposition 5 

is equivalent in this case to the statement that 



T[Y{nY''> -X] 



Q-n-d(Y) 



f5.6l 



1+3 m = -(i(y) + l V-^J 

This is a generalization of the Ward identity (3.5). Note that this identity and 
Proposition 5 hold for arbitrary vertex algebras. 

Since each of the modules Li is automatically a g-module, the Lie algebra gx, 
which was introduced in Sect. 3, acts on the tensor product ®\--yLi. 
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Corollary 1. For any r G J{Li^ . . . , Lp),g G flx and X G ^I^^Li, T{g ■ X) = 0. 

Proof. We have to apply Proposition 5.6 in the case when Y[z) = p[A[z)]^ A G 0. □ 

Corollary 1 means that any Tix^nvariant functional is gx^nvariant. This will allow 
us to prove Theorem 3. 

Proof of Theorem 3. Let Xi{z) and fij{t) be the expansions of (5.1). We assign 
the modules to the points Zi^i = 1, . . . , A^, the modules VF^^(^) to the points 

Wj^j = 1, . . . , m, and the module Wo to the point u of CP^. 

The module Wo is free over the Lie subalgebra of L ® [), generated by a„(m), a G 
A_|_, m < 0; a„(m), a G A_|_, m < 0; hi{m), z = 1, . . . , /, m < 0. In the same way as in 
the proof of Lemma 1, we can show that the space of 7iz^vsr,M^nvariant functionals is 
one-dimensional and is generated by a functional r, whose value on the tensor product 
of the vacuum vectors is equal to 1. It clearly has the property: t(X, = Tjv,m(-^) 
for any X G ®i^iW^A,(^) ®f=\ ^iij(z)-, where v is the generating vector of Wo and Tjv,m 
is the generator of the space J\(z),ijl(z)^ cf. (5.2). 

Recall that we have an embedding p : Vo — ^ Wo of g-modules. Now let Z G Z{q) C 
Vo and consider TN,m{'^-,G,^{-l)v, . . . ,G,^{-l)v ,p{Z)), where to G ®f=iWx^(z) C 
®f=i^\t(z)- This is a linear functional of to. 

We can express this linear functional in two different ways. On the one hand, we 
can represent Z as an element of f/(g_), i.e. in terms of /j[A(m)], A G g,m < 0. By 
Corollary 1 the functional r is invariant with respect to the Lie algebra gz,w,M- So, 
we can use formula (3.5) to "swap" Z. Since Gi^{ — l)v G = 1,... ,m, are 

singular vectors of imaginary weight, A{m) ■ Gi^{ — l)v = for any m > 0. Therefore 
the action of the corresponding elements of the Lie algebra gz,w,M on these vectors is 
equal to 0. So, we obtain, in the same way as in the proof of Proposition 2, 

(5.7) t(cu, G,, (-l)^;, . . . , G,„(-l)^;, ~p{Z)) = [Z{u) ■ , • • • , )](^)- 
On the other hand, we know that p{Z) lies in ttq C Wo. According to (5.6), 

TN,m{i^,G,^{-l)v,. . . ,G,^{-l)v,h,{n)X) = 

rN,m{k{n))TN,m{i^, G,^{-l)v, ... , G,^{-l)v,X). 
By applying this formula several times to p{Z), we obtain 

(5.8) TAr,„(cu, G,, (-l)^;, . . . , G,^{-l)v, p{Z)) = [rN,mMZ)M<.- ■ ■ , <™)(^)- 
From (5.7) and (5.8) we obtain 

Z{u) ■ (P{w\\... ,w^) = [rN,mPz{Z)mw\\... ,w^) 
and Theorem follows. □ 
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Example. For the element S given by (3.4), we have 

Ms) = li:k{-m{-i)-p{-2). 

8 = 1 

Here . . . , /i*, is the basis of [), which is dual to the basis /ii, . . . , /i; with respect to 
the scalar product (•, •) on [) ~ [)* and /? G f) is such that ai{p) = (a^, ai)/2. Therefore 



the eigenvalue of the element S[u) on the vector 



w 



1 1 



w 



is equal to 



Y,h^{u)h*{ 



dup{u) 



1=1 



N 

E 

8 = 1 



A,: 



u 



a,:, 



. ^ u — W,; 
j=i J 



N 



+ E 



Hp) 



E 



For g = sl2 this formula coincides with (2.6). It was written in [20] (cf. formula 
^1.24)). Using this formula we can find the eigenvalues Si[wi ^ . . . , tw^) of the oper- 



ators given by (1.1) on the vector 
eigenvalue is equal to Res 



w 



1 1 



w 



). Since = Res„=^^ S'(m), this 



f5.9l 



Si\W 



1 1 



(u). Thus, we obtain: 



E 



(Aj, Xj) 



E 

s=l 



(A8, ^8 



□ 



Zi 



Lemma 3. The vector cf){w^i , . . . ,1^^) coincides up to a sign with the Bethe vector 
K,... ,w'^) given by (2.8). 

Proof. We will find an explicit formula for 

(5.10) )(^) = TN,m{uJ,G,,V,... ,G,^v). 

and compare it with formula (2.8) for \w\^ ^ . . . , tw^). 

This computation is essentially equivalent to the computation of a bosonic corre- 
lation function from [35]. For X G M* denote by Px the corresponding polynomial 
in Xa^s. Let Px{z) be the local current, which is obtained from Px by replacing Xa 
with a*^{z). In conformal field theory language, 

where Xi G Af^^ , is the bosonic fi^ correlation function 

IN m \ 

(5.11) {\{PxAzi)Y{G.,{w,)) . 

\i=l j=l / 

G 

We will use the generalized Ward identity (5.6) in the case when y = , , where 
Gf)^ is a root generator of the left nilpotent subalgebra n+ of g. 
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This identity reads 

(5.12) = J2 ^iv,™ G-gcu, Gp, i-l)v, ... ,v,... ,GpJ-l)v) + 



c 



+ J2 ^N,m [uj,Gp^{-l)v,. . . ,v,... ,Gp^+pX-'^)v,. . .Gp^{-l)v 



8=1 ^3 \ ■> s / 

where uj G ®^=\^\(z) — ®i^i^A, ^^'^ '^p-y structure constants in the nilpotent 

subalgebra n+ of g. 

By successive use of this identity, we obtain: 

for some co' G (g)fL-^M^^. 

Recall that there is a linear pairing (•,•) : My, X M* — )• C, such that (FiY^X) = 
(y, EiX), z = 1, ...,/, cf. (3.3). Denote by j : Af* — )• C the pairing with the highest 
weight vector of M^: ]{X) = 

When 3^ = a* , the Ward identity (5.6) reads as follows: 

Tjv,m(a^(0)cu,u, . . . ,u) = 0, a G A+. 

Therefore 

TjV,m(cu',U, . . . ,u) = j(cu'). 

If Lo is homogeneous of weight 7, then u' is also homogeneous of weight 7 + Yl'f=i (^ij- 
But = 0, if the weight of cu' is not equal to J2iLi ^i- Therefore ?/'(t/;^\ . . . , w'^ )[(jj) 
0, if the weight of co is not equal to YliLi ~ Yl]Li '^ij- 

Following [35], we obtain by induction the following formula: 

N j{G^...G^ Px,) 

where we used (5.10) and notation from Sect. 3. 
But 

j{G,,...G,^P) = {-irj{E,^...E,,P) 
(cf. Lemma 3.3 from [35]). Therefore we have: 



ijiw\\... ,u;^)(Xi,... ,Xn) 
(5 13) ^ . ..E^jPx^) 



tk, 1 {w-] — W-] ) . . . iw-] — Zi) 
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On the other hand, the pairing of a monomial basis element Fi^ . . . Fi^x^ and 
X G Af* is equal to 

{xy^,F,^...F,,X) =j{F,^...F,,Px). 

Thus, we see that the vector \w^i ■, . . . , w^^) G ®f=\M\^ given by (2.8) defines a linear 
functional on (x)^;^Af^^, which coincides with the functional ( — . . . ,wl^). 
Therefore <p{wW . . . ,wi;^) = {-l)"^\wW. . . ,wi;^). □ 

There is an interesting "analytic" interpretation of Bethe ansatz equations (5.3). 

Proposition 6. The Bethe ansatz equations (5.3) are satisfied if and only if the 
eigenvalue of any Z G Z[g) on the vector (f>[w\^ ^ . . . ,1^^) is non-singular at the 
points Wi^ . . . , ■ 

Proof. Denote by the orthogonal complement to the one-dimensional subspace 
of [) generated by Hi G f). Polynomials in hi[n)^n < 0, form a subspace 7ro,i of 
ttq. Polynomials in h{n),h G ^i~,'n < 0, form a subspace tTq- of ttq. We have: 
7ro(s) = vro,« ® vr^,. 

The space 7ro,i is isomorphic to 7ro(s[2). Denote by Wi C 7ro,i the subspace corre- 
sponding to yV(s[2) C 7ro(s[2). 

Now fix j = 1, . . . , m. The image of Z[g) in ttq is contained in the tensor product 
^o'jj ® [10]. Therefore pz{Z) can be decomposed as J2m XmYm, where G 
'KQ-^,Ym G Wjj- Then we have: 

[rN,mPz]iZ) = J2^N,miXm) ' Tn 

m 

From the definition of the homomorphism rjv,m it is clear that r]\j^ra{Xra) is non- 
singular at M = vOy We will show now that rjv,m(^) is non-singular at u = vo^ for any 
Y G yVij if and only if the Bethe ansatz equation (5.3) is satisfied. 
We know that is the space of differential polynomials in 

P^, = \hi-^f-hi-'^)- 

Since rjv,m is a homomorphism of differential algebras, it is sufficient to check that 
rN,m{Pij) is non-singular at u = if and only if the Bethe ansatz equation is 
satisfied. Applying formula (5.5) we obtain that rjv,m(-P8j) is equal to 

2 V,=i u - z, u-Ws J " V.=i ^ - .=1 u-Ws J ' 

The possible singular terms at u = Wj are 
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which vanishes because ai^[Hi^) = 2, and 

(note that \{Hi.) = 2(A, , )). The latter is non-singular at u = w^, if and 

only if the Bethe ansatz equation (5.3) is satisfied. □ 

Remark 5. Proposition 6 can be considered as a guiding principle for finding eigen- 
values of Gaudin's hamiltonians. This analytic Bethe ansatz has been successfully 
applied to various models of statistical mechanics, cf., e.g., [49]. □ 

In conclusion of this section, let us remark that the Gaudin model can be general- 
ized to an arbitrary Kac-Moody Lie algebra g associated to a symmetrizable Cartan 
matrix. 

Indeed, for such g we can choose the bases {/„} and {/"}, which are orthogonal to 
each other with respect to the invariant scalar product, cf. [6]. Then the operators 
Si given by (1.1) will be well-defined on the A^-fold tensor product of g-modules 
from the category O. For instance, we can take the tensor product of integrable 
representations (/)... (E) ^x^^ or Verma modules (/)... (E) M\j^. 

Moreover, formula (2.8) for the Bethe vector and Bethe ansatz equations make 
perfect sense in this general context. So does formula (5.9). This leads us to the 
following conjecture. 

Conjecture. For an arbitrary symmetrizable Kac-Moody algebra g, the vector 
\w\^ , . . . ,1^^) G (x) . . . (x) given by (2.8) is an eigenvector of the operators 

= 1,... with the eigenvalues Si[wi , . . . ,1^^) given by (5.9), if the Bethe 
ansatz equations (5.3) are satisfied. 

This conjecture can be proved directly using methods of [31], Sect. 7.4-7.8 (note 
that the proof in [31] works for an arbitrary Kac-Moody algebra). One may also 
hope to prove this conjecture by developing theory of Wakimoto modules for general 
Kac-Moody algebras. 



6. Connection with KZ equation. 

We fix A; 7^ —h^. Let io := io{zi, . . . ,zn) be a function with values in V(x). The 
KZ equation is the system of partial differential equations with regular singularities 
on io: 

(6.1) {k + h'')^ = E,uj, z = l,...,iV, 

where S^- is given by formula (1.1). 
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This equation can be considered as the equation on the space -^(x), which was 
defined in Sect. 3 (for more details, cf. [25]). 

Consider the local current S[z)^ corresponding to the element S of the vertex 
operator algebra Vq given by (3.4), 

(6.2) s{z) = i x: : mn^) ■■= + e Luz—'. 

a = l ri£Z 

It is known that the operators Ln generate an action of the Virasoro algebra on any 
representation of g of level k. They have the following commutation relations with 
generators of g: 

(6.3) [Ln,A{m)] = -mA{n + m), A e g. 

Let Cjv be the space with coordinates Zi,... ,Z]\f without the diagonals and 
C'jy be the space C^"*"^ with coordinates Zi, . . . ,Z]\f,t without the diagonals. Denote 
by B and B' the algebras of regular functions on Cjv and C^, respectively. 

Denote by V^^^ the dual space to V^_^j. 

Introduce the S-modules V('^)(z) = V^'^j (^c^, 0jv(z) = gjv^cS, and g'^ = g^c^'. 
The Lie algebra gjv(z) naturally acts on V^^j(z). The Lie algebra g'^ embeds into 
gjv(z) in the same way as the Lie algebra g^ embeds into gjv, cf. Sect. 3. Hence, g'^ 
acts on V^^j(z). The space of invariants of this action, H^^~^{z), is a S-module. By 
Lemma 1, we can identify H^^~^{z) with V^\~^ (E)c ^■ 

Since d/dzi^ z = 1, . . . , A^, are derivations of S, they act on the S-modules gjv(z), 
S;,and Vf^)(z). 

Denote by L^^^i the operator on V^^j(z), which acts as the operator dual to on 
the zth factor and as the identity on all other factors. 

Lemma 4. The operators 

v, = ^-z«;, z = i,... ,iv, 

dz, 

acting on the space V^^j(z), commute with each other and normalize the action of the 
Lie algebra g'^ . 



Proof. We have: 



[V„ V,] = -^xLt + ^L^^: + [L^^L^l*]. 



The first two terms in this commutator vanish, because the operators L^^^ do not 

depend on Zj^ and the last term vanishes, because the operators L^^^ and Z^l* act 
on different factors of V^^j(z). 
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To prove the second statement, consider 

A 



(6.4) 



A G 0. 



We have to show that the commutator of this element with V^, z = 1, . . . , A^, is again 
an element of g^- This will prove the Lemma, because such elements generate g^- 



The action of the element (6.4) on N )^y^{^) is given by 



d 



n-1 



n 



1)1 dz' 



-EE 



m 



s^i m=0 i^J 



\m + l 



(cf. (3.5)). Straightforward computation using (6.3) gives: 

dzi 



—Ai = [M^,L^:},] + 8,,nAi 
Then for the dual operators, acting on V/^^\(z), we obtain 



d 



\J1 r(0* A3*] — A. . „ A3* 

O Zi 



and Lemma follows. □ 



Lemma 4 means that the operators define a fiat connection on the trivial bundle 
over Cjv with the fiber H^x^ ~ ^(A)- Recall that we consider on V^x^ the contragradient 
structure of g-module, with respect to which it is isomorphic to V(x). 

In order to find an explicit formula for this connection, consider for any rj G ^(\) 

the corresponding gz-invariant functional rj G H^^-.. By definition, [L^l^^ ■ rj\{ijj) = 



^(Z« •cu)for anycuGy(x) C Vf;,). 
Using the fact that A(n)(') • e^(cu) 



k + h^ 



The Ward identity (3.5) gives: 
rj{L^li ■ uj) = rj 



k + h' 



'(A)- 

for n > we obtain 

x:/a(-i)«r(o)«-cu. 

a = l 



^^ /'^(0)(0/,(0)0-) J 

j^i a=l 



1 



j^j a = l 



^£ z(/-)(Oz(/,)0-) 



Vk + /i^ 

Thus, the action of on the bundle of invariants is given by 

^ d 
V,: = — 



dz, k + h^ 



30 



BORIS FEIGIN, EDWARD FRENKEL, AND NIKOLAI RESHETIKHIN 



Therefore the flat sections of this bundle are solutions of the KZ equation (6.i). 

Note that this construction can be carried out for arbitrary representations of g 
instead of Va/s. In particular, we can take the contragradient Verma modules Af^^, 
with arbitrary highest weights A^, z = I . . . , A^. Then the flat sections of the bundle 
of invariants will be solutions of the KZ equation with values in (E)iLiMx^. A similar 
analysis of the KZ equation and its generalizations can be found in [26]. 

Schechtman and Varchenko [30, 31] have found integral solutions of the KZ equa- 
tion with values in (E)fLiMx^ and (gj^^M^^. 

We will now derive the first Schechtman- Varchenko solution using Wakimoto mod- 
ules of non-critical level. Essentially this has been done by Awata, Tsuchiya and 
Yamada in [35]. This derivation provides a clear demonstration of how Bethe vectors 
appear in solutions of the KZ equation. 

The construction of the Wakimoto modules from the previous section can be gen- 
eralized to an arbitrary level k [17, 18]. Introduce the Heisenberg Lie algebra [) with 
generators bi[n)^ z = I,... ^l^n ^ Ti^ and the central element 1, with the commutation 
relations 

For X ^ ^* denote by tt^ the Fock representation of [), which is generated by the 
operators bi[n)^n < from the vacuum vector u^, which satisfies 

b,{n)v^ = 0, n > 0, ^^(O)^^ = x{Hi)v^ 

for all z = I, . . . , /, and = v^. 

We define a homomorphism pk from g to the local completion of f^i(r) (x) U{[)) by 
replacing hi{z) with bi{z) in the formulas for the homomorphism p from Sect. 4 and 
shifting by [k -\- h^) the constant q in the formula for p[Fi{z)\. Under the homomor- 
phism pk the central element K G maps to k [17, 18]. Denote the corresponding 
representation of g in Af (x) tt^ by W^^k- 

Let Cp be the space C with coordinates Xi, . . . , without the diagonals and 
be the space C"*"^ with coordinates Xi, . . . , Xp, t without the diagonals. Denote by A 
and A' the algebras of regular functions on Cp and C^, respectively. 

Put W^,fc(x) = 0^1^-;^,^ ®c A rp(x) = Fp (x)c ^ and = {U ® A') ® {U ® 
A'dt) ® ([) (x) A'). The Lie algebra Fp(z) acts on W^^(x); Tix is a Lie subalgebra of 
Fp(x), hence it also acts on Denote by (x) the space of invariants of this 

action. This space is a free module over A, which is generated by the Tix^nvariant 
.4- valued functional {}p on (x)f_;^ VF^^^fc, whose value on Vp := v^^ (/)... (E) v^^ is equal 
to I. 

On the module VF^ we have a natural action of the operators 



(6.5) 



, z — 1, . . . ,p. 
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In the same way as in the proof of Lemma 4 we can show that these operators 
commute with each other and normahze the action of the Lie algebra Tix- Thus, 
these operators act on the space of invariants (x). 

In order to find an exphcit formula for the action of the operators (6.5) on 7?^, we 
have to compute the action of the operator on v\ ^ Tr\ C W\^k- 

It is known that the action of on tta C W\^k coincides with the action of the 
operator 

y y K(nW(-l - n), 

I ragZ r = l 

where ¥[n) are the dual generators to br{—n): 

[¥{n),bs{m)] = nSr^sSn-m^- 
Applying this operator to v\ G tta we obtain 



1 ' 



k + h^ 



r = l 



By Ward identity (5.6), 



1 ^i^^^(,). 



Consider the following system of equations: 
(6.6) ^^^^(X^^_ 

The unique up to a constant factor solution of this system is 

/ = ]^(^^._^^.)(x.x.)/(fc+/^^)_ 

Put Tp(x) = / • This is an Tix^nvariant element of Wj^^(x), which satisfies: 

d 



(6-7) ( - L^^: r^M = 0. 



Now put p = N -\- Xi = Zi^Xi = ^i?^ = I?-- - ^iv+j = ^j^Xiv+j 

—ai^^j = 1, . . . , m. Then we have: 

(6.8) 4,„(z, w)(v^.,„) = Uiz^ - z,f-'^y^'^^') 

i<3 
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where 

hi s<i 

Denote by W\ the subspace of Wx^k^ which is generated from vector v\ by the 
operators a* (0). As a module over the constant subalgebra g of g, W\ is isomorphic 
to M^, cf. Sect. 5. 

The restriction of t^^(z, w) to the subspace 

®L^A, ® G,,{-l)v_^^^ (X) . . . ® G,J-l)v_^^^ 
defines a hnear functional 

This functional vanishes on all homogeneous components except the one of weight 
J2iLi ~ I^jLi ctjj • Therefore it corresponds to a vector 

of weight J2fLi k - EjLi depending on Zi, . . . , zn.Wi, . . . .w^. 

Using the Ward identity (3.5) in the same way as in the proof of Proposition 3, we 
obtain: 

by (6.8). Hence, by Lemma 3, 

(6.9) ^\z, w\\...,w^) = (-1)™ Y{{z, - zi)(^-^^y(^+^''k K , . . . , )• 

Denote by Cm,z the space C™ with coordinates Wi, . . . without all diagonals 
Wj = Wi and all hyperplanes of the form Wj = Zi. The multi- valued function £ defines 
a one-dimensional local system C on the space Cm,z- It is clear that for fixed z the 
value of T^^(z, w) on any vector of the tensor product (E)iLi'^Xt,k ^-a,.,k is a 

product of £ and a rational function on Cm,z- 

Denote by 0(£) = ®i>oO'(£) the twisted de Rham complex of C. It consists of 
differential forms on Cm,z with differential d' acting on a form u by the formula 

d' ■ Lo = d ■ Lo -\- d ■ log L 

It is therefore convenient to represent an element u of the twisted de Rham complex 
by a multivalued differential form I ■ u. The action of d! on u coincides with the 
action oi d on I ■ u. 

If / is a holomorphic function on Cm,z, consider the m-form dwi . . . dwm G 
0™(£). Since it is holomorphic, it defines an element of the mth cohomology group 
H"^{Cm,ZT of the de Rham complex. 
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There is a natural pairing between this group and the mth homology group -ffm(C'm,z, ^*) 
of Cm,z with coefficients in the dual local system. This pairing is given by integration 
over cycles: 

(about different choices of cycles and convergence of the corresponding integrals, cf. 
[51]). This pairing has the following property for any differential form u on Cra,z'- 

(6.10) / d-luj = 0, VA G i?m(C„,z,£*). 

We can now reproduce Schechtman-Varchenko integral formula. 

Theorem 4. Let A he an m-dimensional cycle on Cm,z with coefficients in the dual 
local system C* . The ®f_-^M\^-valued function 



'a 



where \w\^ ^ . . . ,1^^) is the Bethe vector, is a solution of the KZ equation. 

Proof. To avoid confusion, we will denote by A^^'^ or A^'^^'^ the operator, which acts 
as A on the factor Wx^ k or W^a, k of the tensor product ®f-]W\^ k ^-a, fc, and 
as the identity on all other factors. 
From formula (6.7) we obtain: 

^4„(z,w) = zLl)*-4,„(z,w). 

By restricting this equation, we have: 
d 

^4,m(z,w)(cu,G',,(-l)t;_„,^,... ,G,^{-l)v_„^J = 
(6.11) <^^«' 

4,„(z, ■ G,, (-l)^;-„,^ , . . . , G,„(-l)^;_„,^), 

where u G O^^M^^. 

Using the Ward identity, we can rewrite the right hand side of equation (6.11) as 

1 d N 

E E -—-.^N,m{^. w)(/(-)r(^^) • cu, (-l)^;_„,^ , . . . , G.„(-l)^;_„, J 



' a—\ 7=1 ^ J 



-1 d m CO -1 

' a = l s=l ri=0 V ' */ 

r(n)•G,(-l)^;_„,^,,... ,G,„(-l)^;_„,^). 
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To compute the last term of this formula, we have to compute /"(n)-G'i, ( — G 
W_a,^,k for n > 0. 

Lemma 5. For any a = 1, . . . and z = 1, . . . , /, 

r{n) ■ G,{-\)v_^,^ = 0, n>\ 

and there exists such Ya G W-a,,k that 

r{l) ■ G,{-l)v_^, = {k + h'')Ya, r{0) ■ G,{-l)v_^, = {k + h'')L_rYa. 

This Lemma will be proved in the Appendix. It gives for the sth summand of the 
second term of formula (6.12): 

(6.13) 4,„(z,w) (^X:/(-^G,(-l)^;_„,^,... ,G.„(-l)^;_„,^^ , 



where 



According to (6.7), the action of x'^'"-'* on t^^(z, w) coincides with the action of 
d/dws- Therefore we can replace Y^ by 

d Y„. 



dWs Zi — Ws 



Then formula f6.13) can be rewritten as 
d 



dws 



^^4,„(z, w) f E G., (-l)^;-„,^ , . . . , K, . . . , G,„(-l)^;_„,^ 

^« \a = l * / 



For any z and cu, the term in brackets, V, is a product of the multivalued function 
£ and a rational function on C^z- Consider V dwi . . element of 0™(£). 

Then we have 

f d \ ^ 

V dwi . . . dwm = d{V dwi . . . dw^ . . . dwm). 



d 



Ws 



Therefore if A is an m-dimensional cycle on Cm,z with coefficients in £*, we have, 
by formula (6.10): 

( d 



A \ 5 



-V dwi . . . dwyn = 0. 



Ws 



Thus we see that the second term of formula (6.12) will disappear after integration 
over A. 
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But the first term coincides witli tlie action of tlie operator ^i/{k + h^). So, we 
obtain 

J A 



^Jv,m(z,w)(cu,G',,(-i)u_„,^,... ,G',„(-i)u_„,^) dwi...dwr, 



Hence, by (G.ff 
d f 

{k + h)')— / t^^^{z,w){lo,G,^{-1)v_„^^,... ,G,^{-l)v_„^J dwi...dw.. 



'A 

for an arbitrary to G (x)^;^Af^^. 

Tlie last formula can be rewritten as 

{k+K")-^ I -^^z.wW... ,wi^) dwi...dw^ =E, [ *'=(z, , . . . ,wi^) dw^...dw^. 

OZi J A J A 

Therefore we obtain 

{k+h'')^ [ $'=(z,^%... ,w'^) dwi...dw^ = [ ^\z,w\\... ,w'^) dw^...dw^. 

OZi ■' ^ ■' ^ 

Theorem now follows from formula (6.9). □ 

Theorem 4 shows that Bethe vectors enter solutions of the KZ equation. This has 
been observed in [36, 37, 38, 39]. Wakimoto modules provide a tool for constructing 
eigenvectors in Gaudin's model and for solving the KZ equation, thus giving a natural 
interpretation of this phenomenon. 

7. Appendix. 

In this Appendix we give the proof of Lemma 5 and Lemma 2, following [11], 
Appendix B. 

First we prove Lemma 5 in the case g = 5i2- This is a direct computation, using 
the explicit formulas for the homomorphism p (and hence pk) given in Sect. 4. 

We have G{z) = —a{z) = —E[z). From the commutation relations of s[2 we obtain 

[E{n),G{-l)] = 0, [H{n),G{-l)] = 2G{n-l), [F{n),G{-l)] = H{n-l)-k6^,^. 

Since 

E[n)v-a = G[n)v-a = 0,n > 0, H[n)v-a = F[n)v-a = 0,n > 0, 
and H[0)v-a = — 2u_„, we conclude that 

E{n) ■ G{-l)v_a = H{n) ■ G{-l)v_c, = 0, n > 0, 
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and 

F{n) ■ G{-l)v_a = 0, n>l. 

We also have: 

F{1) . Gi-l)v_^, = -{k + 2)v_^,, 

and 

F(0) • G{-l)v_„ = G{-l)F{0)v_„ + H{0)v_„ = 

= (2a(-l)a*(0) - 2a(-l)a*(0) + = -{k + 2)iv_iu_„. 

This proves Lemma 5 in the case g = sl2- 

Recall that in [40], Sect. 5.2, generalized Wakimoto modules over g, associated 
to an arbitrary parabolic subalgebra of g, were defined. Denote by p' the parabolic 
subalgebra of g, which is obtained by adjoining the generator Fi of g to the Borel 
subalgebra of g. Denote by the Lie subalgebra of g generated by Fi^Hi, and Fi. 
We have an orthogonal decomposition of p' into its semi-simple part zi^ and abelian 
part [)^^. 

Denote by L' the Heisenberg Lie algebra with generators a„(n), a* (n), a 7^ G 
Z, and the same commutation relations as in r(g), cf. Sect. 4. Denote by Af' the Fock 
representation of F', which is generated by a vacuum vector satisfying conditions 
(4.1). 

The Wakimoto modules over g of level k associated to p' are realized in the tensor 
product Af' (x)7r^ (x)iv. Here tt^ is the Fock representation of the Heisenberg Lie algebra 

associated to and Z is a representation of 5i2 of level k' such that 
(7.1) k + K" =^-^p^{k' + 2). 

Remark 6. We would like to take this opportunity to correct some errors in [40], Sect. 
5.2. 

First, in the exact sequences on page 179 the module Fndp should be replaced by 
4° 

Second, in formula (4) on page 180 one should take into account the difference in 
normalization of the invariant scalar product (., .) for a Lie algebra and its subalgebra. 
Let Si be the ratio between the scalar product on the Lie algebra g and its Lie 
subalgebra (/W (we use the notation from [40]). Then formula (4) should be rewritten 
as 

k^p ^ + c^(j) = Sikp. 

Since k^ = k -\- Cg, this gives the following condition on levels of representations of 
S-W and g: 

In particular, in our case p = ^\ so 8 = 2/(ai, a^), and we obtain formula (7.1). □ 
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Let us choose as L over SP2, the module Wy.a^^k'^lj £ C. This is a module over 
the Heisenberg algebra generated by aa^{n)^ a*^.[n)^bi[n)^n G Z. It is clear that the 
resulting module over g is isomorphic to the standard Wakimoto module Wy.a^^k^ on 
which the operators G{n) act as — a_„^(n). In other words, we choose such coordinates 
Xa on the big cell U of the flag manifold of g in which Gi = —d/dxa^. In particular, 
we see that the action of the operators Ei[n)^ Hi[n)^ and Fi[n) with n > on the 
subspace Wya^^y C Wy.a^^k is the same as the action of the operators E{n), H{n), 
and F{n) on the s[2-module Wya^^y- 

Hence the action of the operators Ei{n), Hi{n), and Fi[n) on G'i( — I)u_„^ coincides 
with that in the case of sl2- Thus, we obtain: 

E,in) ■ G,i-l)v_^,^ = H,in) ■ G',(-I)t;_„, =0, n > 0, 

E,in) ■ G,i-l)v_^,^ = 0, n > I, 

and 

(7.2) F,(I) • G,(-I)^;_„, = -—^Ak + 

{ rv ■ rv ■ i 



(7.3) F,(0) • G,{-l)v_^, = -{k + /i^)Z_l^;_„,. 

Now consider the operators £'j(n), i7j(n), and Ej[n) with j 7^ z. We first compute 
the finite-dimensional commutation relations of the differential operators Px[Ej]^ Px[Hj]^ 
and PxiF]]} with the vector field Gi (cf. Sect. 4). These commutation relations read 
(cf., e.g., [50], Sect. 3): 

[Px[F,lG,] = 0, [Px[H,lG,] = a,{H,)G,, 

and 

Using these formulas, we can find the commutation relations of Ej[n)^ Hj[n)^ Fj[n)^ 
j 7^ z, with Gi{m), using the Wick theorem. Since Gi[z) = — a„^(2;), there can be 
no "double contractions" , and there is no contribution from the "coboundary term" 
Cjda*^^[z). Therefore this computation amounts to the computation of all "single 
contractions". But those are uniquely defined by the finite-dimensional commutation 
relations above. So we obtain: 

[E,{n), G^z)] = 0, [H,{n), G,{z)] = z-^a,{H,)G,{z), 

and 

[F,{n),G,{z)] = z--a,{H,) : al^{z)G,{z) : . 
These formulas immediately give: 

E,{n) ■ G,{-l)v_„^ = G,{-l)E,{n)v_„^ =0, n > 0, 
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H,{n) ■ G,{-l)v_^,^ = G,{-l)H,{n)v_„^ + a,{H,)G,{n - l)v_„^ = 0, 
F,{n) ■ G,{-l)v_^, = G,{-l)F,{n)v_^, + a,{H,) ^ a:^{r)G,{s)v_^, = 0. 

s-\-r=n — l 

The last two formulas clearly hold for n > 0, and also for n = 0, because i7j(0)u_„^ = 
-a,{Hj)v_a^ and Fj(0)u_„, = -a,(i7j>;^(0)u_„,. 

We can compute the action of any other element of g of the form /"(n), n > 0, on 
G'i( — , using its presentation as a commutator of £'i(m)'s or Fi(m)'s, depending 
on whether belongs to the upper or lower nilpotent subalgebra of g, respectively 
(if G f), then it is a linear combinations of HiS and there is nothing to prove). 
Clearly, in the first case we obtain 0. In the second case, we can realize /"(n) as a 
successive commutator of Ei[n)^n > 0, and Ej[0)^j ^ i. The statement of Lemma 5 
then follows from formulas (7.2), (7.3), and the relation A(0)] = for any A G 0. 

Lemma 2 can be proved along the same lines. We find: 

E,{n) ■ G,{-l)v = H,{n) ■ G,{-l)v = 0, 

and 

F,{n) ■ G,{-l)v = 0,n > 0, F,(0) • G,{-l)v = 8,,,fi^^\H,)v 

for any j = 1, . . . , / and n > 0. Therefore Gi{ — 1) is a singular vector of imaginary 
degree, if and only if ii^^\Hi) = 2{ii^^\ai) I {ai^ai) = 0. 
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